Vulnerability analysis of interdependent infrastructure systems: A methodological framework

Infrastructure systems such as power and water supplies make up the cornerstone of modern society which is essential for the functioning of a society and its economy. They become more and more interconnected and interdependent with the development of scientific technology and social economy. Risk and vulnerability analysis of interdependent infrastructures for security considerations has become an important subject, and some achievements have been made in this area. Since different infrastructure systems have different structural and functional properties, there is no universal all-encompassing ‘silver bullet solution’ to the problem of analyzing the vulnerability associated with interdependent infrastructure systems. So a framework of analysis is required. This paper takes the power and water systems of a major city in China as an example and develops a framework for the analysis of the vulnerability of interdependent infrastructure systems. Four interface design strategies based on distance, betweenness, degree, and clustering coefficient are constructed. Then two types of vulnerability (long-term vulnerability and focused vulnerability) are illustrated and analyzed. Finally, a method for ranking critical components in interdependent infrastructures is given for protection purposes. It is concluded that the framework proposed here is useful for vulnerability analysis of interdependent systems and it will be helpful for the system owners to make better decisions on infrastructure design and protection.     

Spin frustration-related behavior of mixed spin diamond chain

The magnetization of mixed spin-1/2, 1 Ising diamond chain with antiferromagnetic spin exchange interaction has been investigated by using the Monte Carlo simulation. Magnetization behaviors with different antiferromagnetic interactions are simulated under magnetic field at low temperature. The behaviors are studied in two different cases. In one case the absolute value of the interaction between two spin-1/2 atoms is the largest and in the other case it is not the largest. Multiple plateaus of magnetization M, such as plateaus at M=±1, M=0 and M=±0.566 are obtained. The plateaus at M=±0.566 correspond to two metastable states. The M=0 plateau may disappear at some specific values of spin exchange interactions. The spin configurations were checked when the plateaus appear. The dimmer and trimmer states are found.     

GPU‐accelerated direct numerical simulations of decaying compressible turbulence employing a GKM‐based solver

Gas Kinetic Method‐based flow solvers have become popular in recent years owing to their robustness in simulating high Mach number compressible flows. We evaluate the performance of the newly developed analytical gas kinetic method (AGKM) by Xuan et al. in performing direct numerical simulation of canonical compressible turbulent flow on graphical processing unit (GPU)s. We find that for a range of turbulent Mach numbers, AGKM results shows excellent agreement with high order accurate results obtained with traditional Navier–Stokes solvers in terms of key turbulence statistics. Further, AGKM is found to be more efficient as compared with the traditional gas kinetic method for GPU implementation. We present a brief overview of the optimizations performed on NVIDIA K20 GPU and show that GPU optimizations boost the speedup up‐to 40x as compared with single core CPU computations. Hence, AGKM can be used as an efficient method for performing fast and accurate direct numerical simulations of compressible turbulent flows on simple GPU‐based workstations. Copyright © 2016 John Wiley & Sons, Ltd.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

SPH方法研究空气在平板入水中的影响

飞机水上迫降以及船体在高速行驶时,均会发生结构物砰击水面的现象,当结构物底面较为平坦,砰击水面时会捕获空气,在板底形成空气垫。研究空气垫的形成及影响,有助于对结构物进行更准确的载荷分析。该问题的研究涉及到气固液三相的相互作用,在数值模拟中仍是一个挑战,本文使用光滑粒子流体动力学(SPH)方法对该问题进行模拟。首先,与实验结果进行对比;其次,系统地研究了平板形成空气垫以及砰击水面的过程,考查了空气、水和平板的动能变化。对比了在SPH模拟中考虑空气与不考虑空气得到的水压力场和平板所受压力,说明了空气对结果的影响;最后,研究了平板宽度对空气垫形成的影响,从而进一步影响加速度。     

A note on random walks with absorbing barriers and sequential Monte Carlo methods

In this article, we consider importance sampling (IS) and sequential Monte Carlo (SMC) methods in the context of one-dimensional random walks with absorbing barriers. In particular, we develop a very precise variance analysis for several IS and SMC procedures. We take advantage of some explicit spectral formulae available for these models to derive sharp and explicit estimates; this provides stability properties of the associated normalized Feynman–Kac semigroups. Our analysis allows one to compare the variance of SMC and IS techniques for these models. The work in this article is one of the few to consider an in-depth analysis of an SMC method for a particular model-type as well as variance comparison of SMC algorithms.     

Do some chemotaxis-growth models possess Lyapunov functionals?

In spirit of a result by W. Alt from 1980 we give some sufficient criteria that guarantee the existence of Lyapunov functionals for parabolic cross-diffusion models including chemotaxis-growth models with non-diffusive chemotactic signals (resp. with non-diffusive memory).     

Bayesian Estimation for the Order of INAR(q) Mo del

In this paper, we consider the problem of determining the order of INAR(q) model on the basis of the Bayesian estimation theory. The Bayesian es-timator for the order is given with respect to a squared-error loss function. The consistency of the estimator is discussed. The results of a simulation study for the estimation method are presented.     

Synthesis of Depo‐Medrol–chitosan hydrogel as new drug slow‐release appliance and investigation of release kinetics by high‐performance liquid chromatography

The present study deals with preparation and optimization of a novel chitosan hydrogel‐based matrix by suspension cross‐linking method for controlled release of Depo‐Medrol. The controlled release of Depo‐Medrol for effective Rheumatoid arthritis disease has become an imperative field in the drug delivery system. In this context, it was intended to optimize loading circumstances by experimental design and also study the release kinetics of Depo‐Medrol entrapped in the chitosan matrix in order to obtain maximal efficiency for drug loading. The optimum concentrations of chitosan (2.5 g), glutaraldehyde (3.05 μL) and Depo‐Medrol (0.1 mg) were set up to achieve the highest value of drug loaded and the most sustained release from the chitosan matrix. In vitro monitoring of drug release kinetic using high‐performance liquid chromatography showed that 73% of the Depo‐Medrol was released within 120 min, whereas remained drug was released during the next 67 h. High correlation between first‐order and Higuchi's kinetic models indicates a controlled diffusion of Depo‐Medrol through the surrounding media. Moreover, recovery capacity >82% and entrapment efficiency of 58–88% were achieved under optimal conditions. Therefore, the new synthesized Depo Medrol–chitosan is an applicable appliance for arthritis therapy by slow release mechanism. Copyright © 2016 John Wiley & Sons, Ltd.     

Immersed boundary method for unsteady kinetic model equations

Predicting unsteady flows and aerodynamic forces for large displacement motion of microstructures requires transient solution of Boltzmann equation with moving boundaries. For the inclusion of moving complex boundaries for these problems, three immersed boundary method flux formulations (interpolation, relaxation, and interrelaxation) are presented. These formulations are implemented in a 2‐D finite volume method solver for ellipsoidal‐statistical (ES)‐Bhatnagar‐Gross‐Krook (BGK) equations using unstructured meshes. For the verification, a transient analytical solution for free molecular 1‐D flow is derived, and results are compared with the immersed boundary (IB)‐ES‐BGK methods. In 2‐D, methods are verified with the conformal, non‐moving finite volume method, and it is shown that the interrelaxation flux formulation gives an error less than the interpolation and relaxation methods for a given mesh size. Furthermore, formulations applied to a thermally induced flow for a heated beam near a cold substrate show that interrelaxation formulation gives more accurate solution in terms of heat flux. As a 2‐D unsteady application, IB/ES‐BGK methods are used to determine flow properties and damping forces for impulsive motion of microbeam due to high inertial forces. IB/ES‐BGK methods are compared with Navier–Stokes solution at low Knudsen numbers, and it is shown that velocity slip in the transitional rarefied regime reduces the unsteady damping force. Copyright © 2015 John Wiley & Sons, Ltd.